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We put forward a technique to unveil to which symmetry group a nonlinear crystal belongs, making 
use of nonlinear optics with structured light. We consider as example the process of spontaneous 
parametric down-conversion. The crystal, which is illuminated with a special type of Bessel beam, is 
characterized by a nonlinear susceptibility tensor whose structure is dictated by the symmetry group 
of the crystal. The observation of the spatial angular dependence of the lower-frequency generated 
light provides direct information about the symmetry group of the crystal. 
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The spatial arrangements of atoms of substances that can present themselves in a crystalline structure is de¬ 
termined by the specific {point group or symmetry group) category [T] to which the crystal belongs. According 
to Neumann’s principle [2], or principle of symmetry, if a crystal is invariant with respect to a set of symmetry 
transformations, any physical property of the crystal is also invariant under such operations. Thus, the response 
of a crystal to optical, electric or magnetic stimulus can be used to predict general features of its structure, i.e., its 
symmetry group, allowing to unveil its presence in a determined spatial region of interest. 

Crystallographic characterization is usually performed by diffraction techniques using X-rays [3], or material 
waves 0 m that interact with the atoms that are arranged in a certain pattern dictated by the symmetry group. 
By observing the spatial distribution of the photons or particles that come out from the crystal, one can obtain the 
sought-after information. In spite of the high importance of such techniques, sometimes they cannot be used, or its 
use is cumbersome and extremely difficult. Moreover, they are not immune to obstacles that may prevent a unique 
symmetry assignment to a given diffraction pattern. This can be due, e. g., to insufficient number of Bragg peaks 
owing to a finite instrument momentum range or spurious peaks arising from multiple scattering events. In these 
cases, a companion technique might help to identify the appropriate symmetry group of the crystal. 

An alternative to diffraction techniques is the characterization of the crystal symmetry by its effects on the 
nonlinear optical response encoded in the susceptibility tensors [Mi. These tensors are particularly useful because 
of their high sensitivity to both lattice and electronic symmetries. In general, the polarization induced in a material 
when it is illuminated by an optical beam can be written as the Bloembergen expansion [9] 

Pi = + eoX?jlEiE, + e,x%E,E^Ei... ( 1 ) 

where Ei {i = x^y^ z) are the components of the electric field, is the linear susceptibility tensor that determines 

the linear polarization response of the material, and nonlinear susceptibility tensors of different 

ranks responsible for the nonlinear polarization generated in the medium. 

Here we show two main things. Firstly, that it is possible to perform complementary symmetry studies of 
nonlinear materials by using nonlinear optics processes besides second harmonic generation, such as spontaneous 
parametric down-conversion (SPDC). In SPDC, an intense pump beam with frequency uOp interacts with the atoms 

or molecules of a second order nonlinear crystal with nonlinear coefficient In the process, a flow of paired 

photons is generated, the signal and idler, with central frequencies Ug and such that ujg oji = ujp. 

Secondly, that instead of illuminating the nonlinear crystal with paraxial Gaussian beams that propagate along 
a myriad of different directions, as it is usually done in similar cases [7), it is more advantageous to choose as 
illuminating beam a non-paraxial optical beam. The use of a laser source with fast switching of its pointing 
direction might be technically cumbersome, and can severely limit the applicability of the method due to the need 
to control mechanical noise. However, the use of an appropriately designed non-paraxial beam allows to unveil the 
crystal symmetry group in a single shot experiment. 

Under standard conditions cni, most SPDC configurations consist of a Gaussian pump beam that impinges 
normally onto the surface of the nonlinear crystal, which is endowed with a second order nonlinear susceptibility 
Xijk' angle of the crystal is chosen to guarantee the fulfillment of the phase matching conditions, which 

are equivalent to the conservation of energy and momentum of the photons involved in the SPDC process. In most 
cases, the pump beam is linearly polarized. Under these circumstances, the distribution of wave vectors of the 
signal and idler photons, for given frequencies of the resulting photon pairs, is restricted to a single cone (in the 
case of type I SPDC) or to two cones (in the case of type H SPDC) with axes symmetrically arranged relative 
to the pump beam HU. These conical distributions are observed for different crystal symmetries, so they do not 
bear useful information about the point group to which the crystal belongs. This is because in the configuration 
considered, the nonlinear process addresses only certain elements of the nonlinear susceptibility tensor. 

We show that a different configuration can lead to a spatial distribution of the photon pairs that directly reflects 
the symmetry group of the crystal. For the sake of simplicity, we will consider the case of a uniaxial birefringent 
crystal with its optics axis parallel to the normal of its surface. The pump beam that illuminates the crystal is a 
vectorial Bessel mode min], prepared to guarantee a vectorial extraordinary character inside the nonlinear media. 
This can be done by choosing a transverse magnetic (TM) mode with its main propagation direction coincident 
with the axis of the nonlinear crystal. 

A Bessel beam is formed by the superposition of plane waves with wave vectors confined in a cone and with a 
circular cylindrical symmetry on the angular spectrum. In this way, the cylindrical symmetry of the photon pairs 
produced in the SPDC process is directly broken by the intrinsic symmetry of the crystal, which manifests on the 
spatial distribution of the resulting photon pairs. The axicon angle of the Bessel vectorial mode is the parameter 
to be optimized for both the fulfilment of the phase matching conditions and to obtain a clear visibility of the 
characteristic crystallographic pattern built by the photon pairs. 

Note that choosing a Bessel beam as a pump of the nonlinear process is equivalent to observing the crystal 
structure simultaneously for many different angles. In this configuration, SPDC (or any other nonlinear optical 


process) is sensible to all the components of the nonlinear optical susceptibility tensor, including that along the 
main direction of propagation of the pump beam; that is due to the fact that TM modes acquire a significant 
component of its electric field along their main direction of propagation as they depart from the paraxial limit m- 

In this paper we consider vectorial monochromatic beams, i.e., beams with an electric field f (r, t) = 1/2 E^(r) exp 
c.c where uj is the angular frequency of the beam, r designates the spatial location, t is time, is the vectorial 
spatially-varying amplitude of the beam and m is related to the orbital angular momentum of the beam. 

In free space, the electric field of a transverse-magnetic (TM) Bessel beam of order m, which is an exact solution 
of Maxwell’s equations, can be written [12] as a sum of plane-waves with equal amplitude and wavevector confined 
in a cone around the z axis, the main direction of propagation of the beam, with the so called axicon angle cpa^ he., 

E^(r) = ^0 J dpek exp {ik^z -h ip • + imcpp) (2) 

where Sk = k(k • z) — z is the polarization of each wave with wavevector k, the longitudinal wavenumber is 
kz = {oj/c) cos(^a, and the transverse wavevector writes p = (cj/c) sincpa {cos ipp:ksimppf). cpp is the angle 
between the transverse wavevector p and x, and can take any value between 0 and 27r. Notice that even though 
the polarization of each k-wave is perpendicular to k m, the superposition of all k-waves yields a beam with a 
non-zero field component along the main direction of propagation (z). For paraxial beams {cpa small) and m = 0, 
one obtains the so-called radial modes [I2l[l5|. 

In experiments, one never generates such an ideal beam. However, one can generate vectorial Bessel beams close 
to the one described by Eq. 0 by generating superpositions of scalar Bessel modes with the proper topological 
charge and polarization [16]. Under these experimental conditions, the wavevectors k are no longer confined to 
a cone of angle (^a, with transverse wavevector q, but instead the wavevectors spread narrowly around a central 
value, q + Aq, or equivalently, (pa + with |Ag| = \\i\Apa. 

The Hamiltonian of interaction of SPDC is m 

H{t) = Co y" dV J dkp j dkg J dk^ 


Ep(t, r; kp)E^ (t, r; ks)E. (t, r; k^) + h.c. (3) 

where V is the volume of interaction and x^^^ is the nonlinear tensor that characterizes the nonlinear response of 
the material. 

Let us consider as example a type I {eoo) SPDC process in an uniaxial crystal [18]. A configuration that will give 
us information about the symmetry of the crystal is the one where the pump beam propagates inside the nonlinear 
crystal along the optical axis (c = z). The signal and idler waves propagate as ordinary waves. For the case of an 
intense classical extraordinary pump beam that generates ordinary signal and idler photons, the pump beam, and 
the electric field operators for signal and idler photons write m 

E+(t,r;qp) = Eoep(qp) 

E;(t, r; p) = iiV,e,(p) al{K{p), p) ^ 

E-(t,r;q) = iiV,ei(q)al(y(q),q)e*'=^(‘»)^+*‘»-'-, (4) 

where Ng and Ni are normalization factors. 


^*(P) = \J (^) - IpP, kl{q) = /(^) - lql^ 

Cop,os,oi are the ordinary relative permittivities inside the crystal of the pump, signal and idler waves, eep,es,ei are 
the corresponding extraordinary ones, and 




-iUpCop 


kp{kp • z) - z 


es(p)=ksxc, ei(q) = kixc 


( 5 ) 








a|(/c|(p),p) and aj(/c*(p),q) are creation operators for signal and idler photons, with momentum klz + p and 
/cp + q, respectively. 

Since in most situations the nonlinear interaction is weak, we can obtain an accurate quantum description by 
calculating the first-order solution of the Schrodinger equation, i.e., 

|5'(i)) = |vac) - J dt'H{t')\vsic)s\va,c)i (6) 

The interaction Hamiltonian is effectively zero when the classical beams that pump the nonlinear process are zero, 
so that, the time of integration can be extended to t = oc. 

The quantum state of the down-converted photons at the output face of the nonlinear crystal, neglecting for the 
sake of simplicity the contribution from the vacuum term can be written as 

dpdciF{p,q)al{kl,p)al{kl,q)\ya.c) (7) 


where 


^(p,q) = [ep(p + q)]i [es(p)]„ [ei(q)], 

l,m,n 

Afc^(p,q)I/ 


xsmc- 


X exp 


'.kP{p + q) + fc|(p) + fc*(q) j.' 

% - 1j 


(8) 


and Afc^(p,q) = fcP(p + q) - fc|(p) - fc*(q). 

It is important to remind here the crucial role that the tensorial character of plays in the SPDC config- 
uration considered. Most experiments that make use of SPDC, due to the paraxial character of the pump, can 
be described using an effective nonlinear index [m that writes xfff = (ep)j_ (e^)j_ (ei)j_, where are 

the linear polarizations of the pump, idler and signal photons, respectively. The flux of down-converted photons 
depends on the magnitude of the effective index. The vectorial structure of the electromagnetic field defines its 
polarization and, consequently, affects its total angular momentum. Thus, non paraxial pump beams are expected 
to yield interesting results determined by the structure of and its relation to the angular momentum content 
of the down-converted photons m- 

In general, T(p, q) describes an entangled state in the polarization and transverse wavevector degrees of freedom. 
This entanglement takes place both when considering the signal and idler photons as the two subsystems of the 
whole system, and when considering the two degrees of freedom, polarizations and transverse wavevectors, as the 
corresponding subsystems. 

The flux rate i?si(p, q) of detections of a signal photon with transverse wavevector p in coincidence with an idler 
photon with transverse wavevector q [coincidence deteetion)^ or equivalently, a signal photon that propagates inside 
the crystal along the direction {Os^cps) with Og = tan“^ \p\/— and (ps = cos“^Pa,/|p| and an idler 
photon that propagates along the direction [Oi^ipi) with Oi = tan“^ |q|/y^(u;^/c)^eo — |qp and (pi = cos“^ is 


-R(p,q) = l'P’(p,q)P 


(9) 


while the flux rate i?s(p) of signal photons detected with transverse wavevector p [singles deteetion) writes 

-Rs(p) = J dqR{p,q) = J dq|F(p,q)p. (10) 

We illustrate these results, and how they help to determine the symmetry group of the crystal, by showing some 
examples. In Fig. 1, we exemplify the angular dependence of the flux rate of signal photons (singles detections 
as given by Eq. ®) for a variety of nonlinear crystals with different crystallographic symmetries. The numerical 
simulations consider the adequate Sellmeier equations and the crystal birefringent properties. For small values of 
the transverse wave vector of the pump beam (jq^j < 0.01/im“^ for most of the crystals considered here, which 
corresponds to an angle pa ^ 0.0097), the angular spectrum is always formed by a set of concentric rings similar to 
those observed for jp^j < 0.5/im“^ in Fig. 1(a). As jq^j increases the down-converted photons are emitted in a wider 
region in transverse wavevector space, and the symmetry of the crystal becomes more evident. Properties such as the 
radius of the concentric rings, the values of jq^j at which the visibility of the symmetry of the crystal is significant, 
and the maximum value of could be used to determine some characteristics of the linear susceptibility tensor 








and the second- order nonlinear susceptibility tensors. For instance, KDP and CSH 2 ASO 3 belong to the same point 
group, 42 m, but they have slightly different values of e^, Cg and the relevant components of x[mn' manifests in 
general common properties between the angular spectrum of both crystals, but there are still measurable differences 
for small values of due to a higher sensitivity of CSH 2 ASO 3 to the electric field component along the z-axis. 

The coincidence detections, given by Eq. © also provide useful information about the symmetry group of the 
medium under investigation. A direct calculation shows that for jq^j << cj^/c, the detection of a signal photon with 
transverse wavenumber p is accompanied by the detection of an idler photon that is confined to a ring with radius 
|q| ^ |qp|eo/ee around q = — p. As q^ increases, the symmetry of the crystal becomes more visible in the angular 
distribution of idler photons, the flux rate becomes inhomogeneous along the ring with a structure that depends on 
the symmetry of the crystal, and which is compatible with the conditional generation of stationary Bessel photons. 
That is, photons with an angular spectra resulting from superpositions of m and —m Bessel modes. In Fig. 2, this 
effect is illustrated for various types of crystals. 

Summarizing, we have shown that, under the adequate experimental set up, the correlations of the photons 
emitted in an SPDC process, contain crystallographic information that can be accessed by measurements of the 
SPDC angular spectra and the conditional angular spectra. It is expected that a complete characterization of the 
twin photons correlations as a function of the properties of the structured pump beam could be used to obtain 
precise measurements of the first and second order electric susceptibility tensors. Notice that the benefits arising 
from structured pump beams for nonlinear crystallography should not be restricted to the SPDC process. In fact, 
nonlinear harmonic generation crystallography [ 6 ] with Bessel beams could be an alternative to rotational anisotropy 
measurements [7j, since using Bessel modes as a pump is equivalent to simultaneous measurements along different 
directions. 
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FIG. 1. Angular dependence of the flux rate of signal photons detected [singles detection given by Eq. (10)] for type I 
degenerate SPDC (As = = 2Ap). The wavelength of the incident TM Bessel mode is Ap = 407 nm. The angle of the 

cone that characterizes the illuminating pump Bessel beam is sin(y^a = 0.0097. We also consider a spread of directions 
A(/?a ~ 0.0004 rad. The crystal length is taken as L = 1 mm. (a) CdSe: point group 6mm; (b) KDP: point group 42m; (c) 
GaSe: point group 62m; (d) HgGa 2 S 4 : point group 4; (e) HgS: point group 32; (f) LiNbOaipoint group 3m 





















FIG. 2. Flux rate of idler photons detections in coincidence with the signal photons with the transverse wavevector p 
that maximizes Rg. We consider type I degenerate SPDC in a nonlinear crystal of length L = 1mm: (a) GaSe; (b)BBO; 
(c)LiNb 03 and (d) KDP. The incident beam is a Bessel TM mode with wavelength A = 407nm, |qp| = 0.15/im“^ for (a-c) 
and IqpI = 0.2/im“^ for (d). Again, we consider an spread of values of of Aq^ = 1/150/im. 






